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Abstract. We investigate the hydrodynamics of a vari- 
ant of classical Bondi-Hoyle-Lyttleton accretion: a totally 
absorbing sphere moves at various Mach numbers (3 and 
10) relative to a medium, which is taken to be an ideal 
gas having a velocity gradient (of 3% or 20% over one ac- 
cretion radius) perpendicular to the relative motion. We 
examine the influence of the Mach number of the flow and 
the strength of the gradient upon the physical behaviour 
of the flow and the accretion rates of the angular momen- 
tum in particular. The hydrodynamics is modeled by the 
"Piecewise Parabolic Method" (PPM). The resolution in 
the vicinity of the accretor is increased by multiply nesting 
several grids around the sphere. 

Similarly to the 3D models without gradients pub- 
lished previously, models exhibit non-stationary flow pat- 
terns, although the Mach cone remains fairly stable. The 
accretion rates of mass, linear and angular momenta do 
not fluctuate as strongly as published previously for 2D 
models, but similarly to the 2D models, transient disks 
form around the accretor that alternate their direction 
of rotation with time. The average specific angular mo- 
mentum accreted is roughly between 7% and 70% of the 
total angular momentum available in the accretion cylin- 
der and is always smaller than the value of a vortex with 
Kepler velocity around the surface of the accretor. The 
fluctuations of the mass accretion rate in the models with 
small gradients (2%) are similar to the values of the models 
without gradients, while the models with large gradients 
(20%) exhibit larger fluctuations. The mass accretion rate 
is maximal when the specific angular momentum is zero, 
while the specific entropy tends to be smaller when the 
disks are prograde. 

Key words: Accretion, accretion disks - Hydrodynam- 
ics - Instabilities - Shock waves - Methods: numerical - 
Binaries: close 
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1. Introduction 

The simplicity of the classic Bondi-Hoyle-Lyttleton (BHL) 
accretion model makes its use attractive in order to 
roughly estimate accretion rates and drag forces in many 
different astrophysical contexts, ranging from wind-fed X- 
ray binaries (e.g. Anzer & Borner 1996), over supernovae 
(e.g. Chevalier 1996), and galaxies moving through intra- 
cluster gas in a cluster of galaxies (Balsara at al. (1994), 
to the black hole believed to be at the center of our 
Galaxy (Ruffert & Melia 1994; Mirabel et al. 1991). In the 
BHL scenario a totally absorbing sphere of mass M moves 
with velocity relative to a surrounding homogeneous 
medium of density p,^ and sound speed Coo . It has been in- 
vestigated numerically by many workers (e.g. Ruffert 1994 
and 1995, and references therein). Usually, the accretion 
rates of various quantities, like mass, angular momentum, 
etc., including drag forces are of interest as well as the 
properties of the flow, (e.g. distribution of matter and ve- 
locity, stability, etc.). All results pertaining to total accre- 
tion rates are in qualitative agreement (to within factors of 
two, ignoring the instablitites of the flow) with the original 
calculations of Bondi, Hoyle and Lyttleton (e.g. Ruffert & 
Arnett 1994). 

The BHL recipe for accretion in the axisymmetric case 
for pressureless matter is the following. A ring of material 
with radius b (which is identical to the impact parameter) 
far upstream from the accretor and thickness db will be 
focussed gravitationally to a point along the radial accre- 
tion line downstream of the accretor. At this point the 
linear momentum perpendicular to the radial direction is 
assumed to be cancelled. Then, if the remaining energy of 
the matter at this point is not sufficient for escape from 
the potential, this material is assumed to be accreted. The 
largest radius b from which matter is still accreted by this 
procedure turns out to be the so-called Hoyle-Lyttleton 
accretion radius (Hoyle & Lyttleton 1939, 1940a, 1940b, 
1940c; Bondi & Hoyle 1944) 

«» = ^ , a) 
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where G is the gravitational constant. The mass accretion 
rate follows to be 



(2) 



I will refer to the volume upstream of the accretor from 
which matter is accreted as accretion cylinder. 

However, if the assumption of homogeneity of the sur- 
rounding medium is dropped, e.g. by assuming some con- 
stant gradient in the density or the velocity distribution, 
the consequences on the accretion flow remain very un- 
clear. Using the same conceptual procedures, one can cal- 
culate (Dodd & McCrea, 1952; Illarionov & Sunyaev, 1975; 
Shapiro & Lightman, 1976; Wang, 1981) how much angu- 
lar momentum is present in the accretion cylinder for a 
non-axisymmetric flow which has a gradient in its density 
or velocity perpendicular to the mean velocity direction. 
Then, assuming that the angular momentum will be ac- 
creted together with the mass, it is only a small step to 
conclude that the amount of angular momentum accreted 
is equal to (or at least is a large fraction of) the angular 
momentum present in the accretion cylinder. Note, that 
if the velocity is a function of position, then by virtue of 
Eq. (Q) also the accretion radius varies in space. Thus the 
cross section of the accretion cylinder (perpendicular to 
the axis) is not circular. 

However, the reasoning of BHL calls for a cancelling 
of linear momentum perpendicular to the radial accre- 
tion line before matter is accreted. Together with this 
linear momentum also angular momentum is cancelled 
and so the matter accreted has zero angular momen- 
tum by construction! This point was first discussed by 
Davies & Pringle (1980), who were able to construct two- 
dimensional flows with small non-vanishing gradients for 
which the accreted angular momentum was exactly zero, 
by placing the accretion line appropriately. Thus, follow- 
ing these analytic investigations two opposing views are 
voiced about how much angular momentum can be ac- 
creted: either a large or a very small fraction of what is 
present in the accretion cylinder. Numerical simulations 
thus are called for to help solve the problem. 

In this paper I would like to compare the accretion 
rates of several quantities (especially angular momen- 
tum) of numerically modeled accretion flows with gradi- 
ents to the previous results of accretion without gradi- 
ents (e.g. Ruffert 1994). One has to change some of the 
parameters of the flow (Mach number, size of the accre- 
tor) in order to get a good overview of which features are 
generic and which specific to that combination of param- 
eters. Although several investigations of two-dimensional 
flows with velocity gradients exist (Anzer et al. 1987; Fryx- 
ell & Taam 1988; Taam & Fryxell 1989; Ho et al. 1989), 
i/iree-dimensional simulations are scarse due to their in- 
herently high computational load. Livio et al. (1986) first 
attempted a three-dimensional model including gradients, 
but due to their low numerical resolution the results were 
only tentative. Also in the models of Ishii et al. (1993) was 



the accretor only coarsly resolved, while the results of Bof- 
fin (1991) and Sawada et al. (1989) are only indicative, be- 
cause due to the numerical procedure the flows remained 
stable (too few SPH particles in Boffin 1991 and local time 
stepping in Sawada et al. 1989 which is appropriate only 
for stationary flows). A simulation that was numerically 
better resolved was performed later by Ruffert & Anzer 
(1995), but since only one model was presented, the results 
cannot be taken as conclusive either. I intend to remedy 
these shortcomings in the present paper. 

In section || I give only a short summary of the numer- 
ical procedure used. Sections |J to^ present the results, 
which I analyze and interpret in Sect.|8[ Section ^ sum- 
marizes the implications of this work. 

2. Numerical Procedure and Initial Conditions 

Since the numerical procedures and initial conditions are 
mostly identical to what has already been described and 
used in previous papers (e.g. Ruffert, 1996; Ruffert & 
Anzer, 1994) I will refrain from repeating every detail, 
but only give a brief summary. 

2.1. Numerical Procedure 

The distribution of matter is discretised on multiply 
nested equidistant Cartesian grids (e.g. Berger & Colella, 
1989) with zone size S and is evolved using the "Piece- 
wise Parabolic Method" (PPM) of Colella & Woodward 
(1984). The equation of state is that of a perfect gas with 
a specific heat ratio of 7 = 5/3 or 7 = 4/3 (see Table [[]). 
The model of the maximally accreting, vacuum sphere in a 
softened gravitational potential is summarized in Ruffert 
& Arnett (1994) and Ruffert & Anzer (1995). 

A gravitating, totally absorbing "sphere" moves rela- 
tive to a medium that far upstream has a distribution of 
density and velocity given by 



Vq [ 1 + — tanh 



2e v 



Ra 



0, v z 



(3) 



0,(4) 



with the redefined accretion radius 

2GM 
Ra = — o - ■ 



(5) 



In this paper I only investigate models with gradients 
of the velocity distribution; the values of e v can be found 
in Table 0. Thus for all models I set e p = 0. Additionally, if 
only a density gradient is introduced without varying some 
other thermodynamic variable (e.g. temperature, entropy, 
etc.) at the same rate, pressure will not be in equilibrium 
(cf. e.g. Ho et al. 1989), so an additional thermodynamic 
variable should be varied, which complicates matters. 
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Table 1. Parameters and some computed quantities for all models. Mac is the Mach number of the unperturbed flow, e v the 
parameter specifying the strength of the gradient, 7 the ratio of specific heats, _R* the radius of the accretor, g the number of 
grid nesting depth levels, 8 the size of one zone on the finest grid, e the softening parameter (zones) for the potential of the 

accretor (see Ruffert, 1994), if the total time of the run (units: Ra/coo), M the integral average of the mass accretion rate, S 

one standard deviation around the mean M of the mass accretion rate fluctuations, M the maximum mass accretion rate, Mbh 
is defined in Eq. (3) of Ruffert & Arnett (1994), i x , l y , l z , are the averages of specific angular momentum components together 
with their respective standard deviations <r x , a y , <r z , s is the entropy (Eq. (4) in Ruffert & Arnett 1994), the number TV of zones 
per grid dimension is 32, and the size of the largest grid is L — 32_Ra (except for model RL for which it is L = 128-Ra). 
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3 


-0.03 


5/3 


0.02 


10 


1/512 
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4.82 


0.72 


± 


0.04 


0.78 


0.00 ± 


0.01 


0.00 


± 


0.02 


+0.20 ± 0.04 


2.1 


IS 


3 


-0.03 


5/3 


0.02 


9 


1/256 


3 


13.9 


0.53 


± 


0.09 


0.80 


0.00 ± 


0.13 


0.05 


± 


0.17 


+0.12 ± 0.24 


2.2 


IM 


3 


-0.03 


5/3 


0.10 


7 


1/64 


4 


26.6 


0.79 


± 


0.06 


0.90 


-0.01 ± 


0.11 


0.00 


± 


0.13 


+0.68 ± 0.33 


2.2 


IM* 


3 


-0.03 


5/3 


0.10 


7 


1/64 


4 


8.37 


0.82 


± 


0.08 


0.91 


+0.02 ± 


0.06 


-0.08 
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0.14 


+0.69 ± 0.45 


2.2 


JS 


10 


-0.03 


5/3 


0.02 


9 


1/256 


3 


2.93 


0.45 
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0.09 


0.68 


-0.04 ± 


0.17 


-0.08 
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0.35 


+0.18 ± 0.28 


5.3 


JM 


10 


-0.03 


5/3 


0.10 


7 


1/64 


4 


10.3 


0.72 


± 


0.05 


0.79 


+0.01 ± 


0.06 


-0.07 
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0.33 


+0.26 ± 0.35 


5.2 


KS 
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5/3 


0.02 
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1/256 
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6.89 


0.54 
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0.09 


0.78 


0.00 ± 


0.02 


-0.01 
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0.03 


+0.07 ± 0.02 


1.8 


KM 
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-0.20 


5/3 


0.10 
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1/64 


4 


20.3 


0.95 
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0.19 


1.30 


0.00 ± 


0.02 


0.00 
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0.02 


+0.26 ± 0.06 


1.6 


LS 


10 


-0.20 


5/3 


0.02 
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1/256 
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1.94 


0.35 
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0.08 


0.53 


+0.02 ± 


0.03 


+0.01 
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0.05 


+0.09 ± 0.03 


5.2 


LM 


10 


-0.20 


5/3 


0.10 
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1/64 
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8.54 


0.72 


± 


0.17 


1.11 


0.00 ± 


0.04 


0.00 


± 


0.05 


+0.25 ± 0.09 


4.7 


RL 


0.6 


-0.20 


5/3 


1.00 


6 


1/8 


5 


63.1 


36.3 


± 


0.14 


36.4 


0.00 ± 


0.00 


0.00 


± 


0.00 


-0.49 ± 0.04 


0.18 


ST 


3 


-0.03 


4/3 


0.02 


10 


1/512 


8 


4.60 


1.01 


± 


0.09 


1.29 


0.00 ± 


0.05 


+0.01 


± 


0.07 


+0.51 ± 0.23 


4.4 


SS 


3 


-0.03 


4/3 


0.02 


9 


1/256 


3 


9.24 


1.01 


± 


0.12 


1.46 


+0.02 ± 


0.15 


+0.03 


± 


0.25 


+0.36 ± 0.40 


5.2 



The function "tanh" is introduced in Eq. (Q) to serve as 
a cutoff at large distances y for large gradients e v : In some 
models I imposed a gradient of e v — 0.2 which at distances 
beyond 5Ra would produce negative velocities if a linear 
distribution were used. In the limit of small y <C Ra, 
Eq. (U) transforms to a shape very similar to Eq. (^) with 
p replaced by v (as in Ruffert & Anzer 1995, Eq. (2)). 
The relative velocity vq is varied in different models: I 
perform simulations with Mach numbers Moo = ^o/coo 
of 3.0 and 10. In the reference frame of the accretor the 
surrounding matter flows in +x-direction. Our units are 
(1) the sound speed Coo as velocity unit; (2) the accretion 
radius (Eq. (||)) as unit of length, and (3) poo = po as 
density unit. Thus the unit of time is Ra/coo- 

When using the density and velocity distributions 
(Eqs. (||) and (Q)) to calculate the mass accretion rate, 
assuming that all mass within the accretion cylinder is 
accreted, one obtains to lowest order in e v 

M = irRlpoVo , (6) 

an equation very similar to Eq. (|2|). Further assuming 
that all angular momentum within the deformed accretion 
cylinder is accreted too, the specific angular momentum of 
the accreted matter follows to be (Ruffert & Anzer 1994; 
Shapiro & Lightman 1976; again to lowest order in e v ) 

jz = j (6e v - £ P ) v Ra ■ (7) 



Note the different signs with which the two e enter this 
equation. The density gradient e p acts in the "expected" 
direction: if the density is higher on the positive side of the 
y-axis, then the vortex formed around the accretor is in 
the counter-clockwise direction, i.e. the angular momen- 
tum is negative. Contrary to this, if the velocity is larger 
on the positive y-side, then the shortened accretion radius 
on this side reduces the cross-section for the higher specific 
angular momentum to such an extent that the rotational 
direction of the vortex is reversed: the angular momentum 
is positive. 

The values obtained from the numerical simulations of 
the specific angular momentum should be compared to the 
values that follow from this Eq. (Q) to conclude which of 
the above mentioned views — low or high specific angular 
momentum of the accreted material — is most probably 
correct. Ruffert &+^.nzer (1995) find that the numerical 
value is 0.72 timestl the analytical estimate, which would 
indicate a large value for the accreted specific angular mo- 
mentum. Although the simulations yield the values of all 
three components of the angular momentum, the interest- 
ing component is the one pointing in z-direction. Thus I 
will implicitely assume j z when discussing properties like 
fluctuations, magnitudes, etc. that the simulations pro- 
duce. From the symmetry of the boundary conditions the 



this value is the rms, while in Table hi we list the averages 
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Fig. 1. The coefficient for the mass accretion rate / m (dotted), 
defined by Eq. (|^), and for the specific angular momentum /j 
(solid), defined by Eq. (^|), as a function of e v - The two thin 
curves that end at around -0.8 show the values of / for a simple 
linear relation between w xoo and e v , while the bold curves apply 
to a relation including the "tanh"-term as given in Eq. (||). The 
two vertical straight lines indicate the gradients that were used 
in the numerical models (cf. Table |l|), e v = 0.03 and e v = 0.2. 



average of the x and y components of the angular mo- 
mentum should be zero, although their fluctuations can 
be quite large. Numerically I obtain the mass and angular 
momentum accretion rates as a function of time, M(t) and 
J(t). From these functions I calculate the instantaneous 
specific angular momentum j(t) = J(t)/M(t) as function 
of time, which I will plot and from which I calculate the 
temporal mean I listed in Table § 

One can numerically approximate the integrals 
(U. Anzer, personal communication) of the mass flux and 
angular momentum over the deformed cross section of the 
accretion cylinder, to obtain the coefficients / in the rela- 
tions Eq. (§) and Eq. (g): 



M = / m (ev)7ri?|poWo 



.h = fj(e v )e v v R A 



(8) 



(9) 



Here, I will only consider the effect of a velocity gradi- 
ent. The unitless functions / are a function of e v and the 
functional relation of v X oo[£v], i.e. whether v xoo depends 
purely linearly on e v or as in Eq. (|J) via the "tanh"-term. 
Figure [l] shows the values of the functions / for the mass 
and specific angular momentum and for both the linear 
and "tanh" case. In the linear case there is no solution 
for e v 0.15 (U. Anzer, personal communication) so the 
curves end at that point. Since / m w 1 and is practically 



constant for e v ^ 0.1, the Eq. (j6[) is a good approximation 
in this range. If the prescription is correct that everything 
in the accretion cylinder is accreted, we expect to see an 
increase of the mass accretion rate by a factor of roughly 
1.8 in the models with a fairly large gradient of e v = 0.2 
(cf. Table 0). The same trends apply to the specific an- 
gular momentum: its coefficient /j remains relatively con- 
stant /j « 1.5 in the range s v ^ 0.1, and becomes a factor 
of 2 larger for e v s» 0.2. In the case including "tanh", the 
cocfficent /j decreases again for e v <; 0.25 because the 
gradient is so steep that the "tanh" -cutoff acts at very 
short distances. So the short lever arm that enters into 
the angular momentum wins. 

2.2. Models 

The combination of parameters that I varied, together 
with some results are summarised in Table [j] The first 
letter in the model designation indicates the Mach num- 
ber and the strength of the gradient: I, J, and S have 
e v = 0.03, while K, L, and R have e v = 0.2. The second 
letter specifies the size of the accretor: M (medium) and 
S (small) stand for accretor radii of 0.1 Ra, and 0.02 Ra, 
respectively. I basically simulated models with all possible 
combinations of two relative wind flow speeds (Mach num- 
bers of 3 and 10), two gradient strengths (3% and 20%) 
and two different accretor sizes (0.02 and 0.1 accretion 
radii), all with an adiabtic index of 5/3. The exeptional 
models are ST and SS — in which I used an adiabatic 
index of 4/3 — and model RL which has a very large ac- 
cretor radius and a very slow relative flow velocity. The 
grids are nested to a depth g such that the radius of the 
accretor R+ spans several zones on the finest grid and the 
softening parameter e is then chosen to be a few zones 
less than the number of zones that the accretor spans. 
Model IT and ST are physically identical to models IS 
and SS, respectively. However, models IT and ST are nu- 
merically better resolved, because they are nested one grid 
level finer. Model IM is identical to the model presented in 
Ruffert & Anzer (1995), i.e. the velocity gradient is cho- 
sen as specified in Eq. (2) of Ruffert & Anzer (1995) with- 
out the "tanh" -term of Eq. (|J). This term is included in 
model IM*, however, since the results were nearly indis- 
tinguishable between models IM and IM*, model IM* was 
evolved for only roughly one third of the time of model IM. 

As far as computer resources permitted, I aimed at 
evolving the models for at least as long as it takes the flow 
to move from the boundary to the position of the accretor 
which is at the center (crossing time scale). This time is 
given by L/2Moo and ranges from about 1 to about 10 
time units. The actual time if that the model is run can 
be found in Table |l|. 

The velocity distribution following the "tanh"- 
prescription of Eq. (|J) has an inflection point and thus 
is Kelvin-Hclmholz unstable with an amplification time 
constant of roughly t ps (5Ra) / {SvMcoo) (Drazin 1981). 
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Fig. 2. Contour plots showing snapshots of the density to- 
gether with the flow pattern at large distances from the accre- 
tor. The contour lines are spaced logarithmically in intervals 
of 0.1 dex. The bold contour levels are sometimes labeled with 
their respective values (0.01 and 0.4). The dashed contour de- 
limits supersonic from subsonic regions. The time of the snap- 
shot together with the velocity scale is given in the legend in 
the upper right hand corner of each panel. 
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During the time it takes matter to move from the bound- 
ary to the accretor, t « (16i?A)/(A / ic 00 ), random pertur- 
bations can grow by about t/r 3e v which is smaller 
than unity even for the large gradients (e v = 0.2) used in 
the simulations and listed in Table [I]. 

The calculations are performed on a Cray-YMP 4/64 
and a Cray J90 8/512. They need about 12-16 MWords of 
main memory and take approximately 40 CPU-hours per 
simulated time unit (for the 6 = 1/64 models and Mach 
10; the 8 — 1/256 models take four times as long, etc.; 6 
is the size of a zone on the finest grid, see Table [[]). 

3. Shape of the shock cone 

Figure ^ shows the density distribution of five models to- 
wards the end of the simulations, emphasizing the distri- 
bution of matter on scales between one and ten accretion 
radii. The left Figs. || a, c and e, show models with a small 
gradient in velocity (3%), while the right Figs, b and d 
have large gradients (20%). Similarly to 3D models with- 
out gradients (see Ruffert 1996, and references therein), 
these new models do not exhibit the "flip-flop" flow vis- 
ible in previous 2D simulations. The shape of the shock 
cones shown in Fig. || is fairly constant in time and re- 
mains roughly conical, contrary to the 2D flows, whose 
cones shifted strongly from side to side. One notices the 
following features when inspecting Fig. |^. The mass is dis- 
tributed in a hollow shock cone (as has been reported pre- 
viously) for the models with a small gradient, i.e. the den- 
sity is maximal just behind the shock, while downstream 
from the accretor, the density is minimal along the axis. 
The asymmetry of the velocities in the incoming flow re- 
flects itself in higher density maxima along the cone on the 
side of the lower velocities. This density asymmetry is so 
pronounced in the models with strong gradients (Figs. ||) 
that the "hollow cone" shape can be recognized only with 
difficulty. The line of minimum density is very irregular 
and is shifted from the y = 0-axis by several accretion 
radii. Already upstream of the shock a higher density is 
indicated by the contour of value 0.01 being detached from 
the shock on the positive y-axis side, while it is very close 
to the shock on the negative y-axis side. This density dif- 
ference is easily explicable: on the side of smaller velocities 
gravity can act relatively more strongly to divert the flow. 
Thus the effective local accretion radius is larger on the 
side of smaller velocities and so a larger volume of matter 
can be gravitationally focused by the accretor. One also 
notices that the side of the cone with smaller densities is 
more irregularly shaped than the high-density side. The 
cavities and lumps produced by the fluctuating flow close 
to the accretor (at distances closer than roughly one ac- 
cretion radius) can propagate more easily downstream on 
the side of the cone with lower densities. Since the ve- 
locity enters the accretion radius Eq. |J) via a square, one 
might wonder, whether the velocity is so small, that the 
local accretion radius is comparable to the distance of the 



accretor to the boundary of the computational box which 
is approximately at 16i?A- This is not the case, since in- 
serting Vrain = v /2 (from Eq. (0)) into Eq. ||) one obtains 
4i?A, which is a factor of 4 smaller than the distance to 
the computational boundary. 

4. Results of models with 3% velocity gradient 

4-.1. Moderately supersonic accretors, Mach 3 

The model denoted by IM in this work is identical to the 
simulation presented in Ruffert & Anzer (1995). Since for 
the parameters that I wanted to investigate in this paper 
a slightly different prescription of the velocity gradient 
was necessary (cf. Sect. |2.2[ ), a comparison is called for to 
check how large the influence of the "tanh"-term is for 
small gradients e v — 0.03. Models IM and IM* are iden- 
tical except for the "tanh"-term. Fig. ||a shows a contour 
plot of the density in one plane for model IM, Fig. |^ shows 
the accretion rates of several quantities for model IM, and 
Fig. H displays the same quantities for model IM*. One 
notices that both the mass accretion rate as well as the 
specific angular momentum accreted are similar in both 
models for the time over which both have been calcu- 
lated. I stopped the simulation of model IM* at t » 8.2 
because it looked so alike to model IM. Thus I conclude 
that for small gradients e v = 0.03 the difference between 
the "tanh" -prescription and the simple linear dependence 
is indeed negligible. 

Models IS and IT differed only in their numerical res- 
olution: model IT was simulated with one grid level finer, 
while the accretor size was very small (R — 0.02i?A) com- 
pared to models IM and IM* (whose radii were 5 times 
larger). The density distribution of models IS and IT can 
be found in Fig. ||, while their accretion rates are shown in 
Fig. ||. Because model IT had one level of refinement more, 
the computational cost was larger, so the time the simula- 
tion was run is roughly half of the time of model IS. Until 
(«5 both models show the same features: the mass ac- 
cretion rate rises to roughly 6 units and starts fluctuating 
at about i«4 time units. The specific angular momen- 
tum rises continuously until t ~ 1 then fluctuates around 
a value of roughly -0.03 while increasing the amplitudes 
of the fluctuations. Thus the use of 9 nested grids seems 
sufficient, since the calculation with 10 grids did not show 
any qualitative differences. 

It is clear that model IT has not been evolved for long 
enough to obtain meaningful averages, since the fluctua- 
tion of the mass accretion rate has hardly begun when the 
simulation is stopped (Fig. [|). Although model IS was cal- 
culated for a longer time, the continuously decreasing mass 
accretion rate of model IS (Fig. ||) indicates that a station- 
ary state has not yet been reached und so the simulation of 
this model should have been continued even further. How- 
ever, the high computational cost precluded this. The fluc- 
tuations of the specific angular momentum, too, seem to 
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Fig. 3. Contour plots showing snapshots of the density to- 
gether with the flow pattern in a plane containing the center of 
the accretor for all models with a velocity gradient of 3%. The 
contour lines are spaced logarithmically in intervals of 0.1 dex. 
The bold contour levels are sometimes labeled with their re- 
spective values (0.01 and 1.0). The dashed contour delimits 
supersonic from subsonic regions. The time of the snapshot 
together with the velocity scale is given in the legend in the 
upper right hand corner of each panel. 



increase with time corroborating the the statement. Thus 
the average quantities obtained from this model IS cannot 
be very exact. 

Two uncertainties that could not be resolved by the 
single model IM presented in Ruffert & Anzer (1995), can 
now be answered. The first one pertained to the fact that 
the specific angular momentum (visible in the top right 
panel of Fig. |J) seemed to reach but not cross the zero- line. 
This seems to be a coincidence of the initial and boundary 
conditions, since in model IS (visible in the middle right 
panel of Fig. |j) the fluctuations are indeed large enough 
to change the sign of the specific angular momentum ac- 



creted. We will return to this point in Sect. 3.2. The sec 



ond uncertainty was whether in the generic case the values 
of the specific angular momentum attained and exceeded 
the analytically estimated ones given byEq. (|). Model IS 
clearly does not attain these values by a large margin, 
roughly a factor of 3 — the analytic value for model IS is 
the same as for model IM, since the accretor radius does 
not enter into Eq. (]?]). The smaller accretor radius seems 
to allow only smaller values of angular momentum to be 
accreted: if the lever arm (which is the radius of the ac- 
cretor) is smaller the velocities have to be an appropriate 
amount larger (a factor of 5, roughly) to compensate. This 
is obviously not the case: the arrows in the left panels of 
Fig. H close to the surface of the accretor have roughly 
the same length. The smaller accretor sizes also have the 
effect that the time scale of the fluctuations of model IS 
are shorter than the ones of model IM (compare the right 
panels of Fig. ||) . 

The corresponding plot to the top left panel in Fig. |I] 
of model IM for the axisymmetric case can be found in the 
top left panel of Fig. 16 in Ruffert & Arnett (1994). For 
model IS the closest would be the top left panel of Fig. 22 
in Ruffert & Arnett (1994). One can see, that while the 



amplitude of the fluctuations of the the z-component of 
the accreted angular momentum is comparable, the av- 
erage of this component of the models with velocity gra- 
dients is clearly non-zero. Contrary to this, the x- and y- 
components fluctuate more strongly in the models without 
gradients, but in all models their temporal average is close 
to zero (see e.g. Table [j]). The run, average and fluctua- 
tions of the mass accretion rate is similar in all models. 

Due to the non-axisymmetric upstream boundary con- 
ditions it is not surprising that the shape of the bow shock 
is not symmetric either. There is an indication of this fact 
in the panels shown in Fig. |^ for the density distribution 
close to the accretor, but it is very prominent when in- 
specting the shock cone position further away from the 
accretor (see e.g. Fig. 1 in Ruffert & Anzer 1995). The 
temporal evolution shows the usual kinks and deforma- 
tions of the shock cone that were described in the previous 
papers (Ruffert 1996, and references therein). 

4-- 2. Highly supersonic accretors, Mach 10 

The right panels of Fig. ||| display the density contours 
of models JM and JS, which are equivalent to mod- 
els IM and IS, except for the different flow speed upsteam: 
Mach 10 for the J-models contrary to Mach 3 for the I- 
models. The corresponding accretion rates can be found 
in Fig. g. The highly supersonic models JM and JS, too, 
do not converge to a quiescent steady state but show an 
unstable fluctuating flow pattern. 

In the same way as model IM, also model JM meets 
and exceeds the analytically given value (Eq. (§)) of the 
specific angular momentum (top right panel Fig. |6J, hori- 
zontal line), however only in rare, short bursts. Thus, on 
average, the fraction of the specific angular momentum to 
the analytic value is smaller than the value of the fraction 
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Fig. 4. The accretion rates of several quantities are plotted as a function of time for the moderately supersonic (A^oo=3) 
models IM, IS and IT with a velocity gradient of 3%. The left panels contain the mass and angular momentum accretion rates, 
the right panels the specific angular momentum of the matter that is accreted. In the left panels, the straight horizontal lines 
show the analytical mass accretion rates: dotted is the Hoyle-Lyttleton rate (Eq. 1 in Ruffert 1994a), solid is the Bondi-Hoyle 
approximation formula (Eq. 3 in Ruffert 1994a) and half that value. The upper solid bold curve represents the numerically 
calculated mass accretion rate. The lower three curves of the left panels trace the x (dotted), y (thin solid) and z (bold solid) 
component of the angular momentum accretion rate. The same components apply to the right panels. The horizontal line in 
the right Panel of model IM shows the specific angular momentum value as given by Eq. (Q) . It is outside the range of the plot 
for models IS and IT. 
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Fig. 5. The accretion rates of several quantities are plotted as a function of time for model IM*. The left panel contains the 
mass and angular momentum accretion rates, the right panel the specific angular momentum of the matter that is accreted. In 
the left panel, the straight horizontal lines show the analytical mass accretion rates: dotted is the Hoyle-Lyttleton rate (Eq. 1 
in Ruffert 1994a), solid is the Bondi-Hoyle approximation formula (Eq. 3 in Ruffert 1994a) and half that value. The upper 
solid bold curve represents the numerically calculated mass accretion rate. The lower three curves of the left panels trace the 
x (dotted), y (thin solid) and z (bold solid) component of the angular momentum accretion rate. The same components apply 
to the right panel. The horizontal line in the right Panel shows the specific angular momentum value as given by Eq. (Q). 



of model IM (cf. Table [l] colum l z ). However, similarly to 
the difference between models IM and IS, the specific an- 
gular momentum of model JS is smaller still. In fact, when 
looking at the bottom right panel of Fig. |^, the curve for 
the j z component overlaps very strongly with the curves 
of the other components. There is, however, still clearly a 
systematic shift of j z not visible in j x and j y which fluc- 
tuate around zero. 

The left two panels in Fig. |] (models JM and JS) can 
be compared to the two top panels of Fig. 9 in Ruffert 
(1994b) which show the results of the equivalent models 
without gradients (FM and FS, respectively). The mag- 
nitude and fluctuation amplitude of the mass accretion 
rate are similar which is also the case for the angular mo- 
mentum accretion rates (note different scale of y-axis in 
figures). Since the mass accretion rate of model JS seems 
to be steadily declining (although also fluctuating, bottom 
left panel in Fig. |6|) , this model at the time we stopped the 
simulation is probably still influenced by transients gen- 
erated through the initial conditions. Thus the average 
values of e.g. the specific angular momentum, have to be 
used with caution. 

5. Results of models with 20% velocity gradient 

Models KM, KS, LM and LS break the axisymmetry very 
strongly since they have as boundary condition a 20% ve- 
locity gradient over one accretion radius. This asymmetry 
induces a strong clockwise vortex around the accretor as 
can be seen in Fig. (?]. The sign of the velocity gradient 



has been chosen in such a way that upstream of the accre- 
tor the higher velocities are on the negative y-side of the 
xy-plane (lower half of the contour plots). However, the 
numeric simulations confirm the sign of the analytically 
estimated accreted angular momentum (Eq. (0)): the vec- 
tor points into the plane of the plot, which corresponds to 
a clockwise rotation in these contour plots. A note of cau- 
tion is necessary, however: recall that Eq. was derived 
only to lowest order in e v and e p , thus assuming that e v 
and e p are small compared to unity. Obviously this is ques- 
tionable with the choice e v — 0.2 for the models presented 
in this section. It is still interesting to see by how much 
the analytic estimates deviate from the numerical models 
is this extreme case. 

5.1. Moderately supersonic accretors, Mach 3 

The accretion rates of several quantities for models KS 
and KM can be found in Fig. |[ Contrary to model IM, 
model KM does not exhibit such strong fluctuations, thus 
it never even comes close (by a factor 3) to the analyt- 
ically predicted value (-0.9) of Eq. (^), nor does it ever 
come close to the zero line. The large scale motion of the 
vortex around the accretor is fairly stable in time, which 
explains why the fluctuations relative to the mean of the 
specific angular momentum are smaller in model KM than 
in model IM. 

The fluctuations of the mass accretion rate (left pan- 
els of Fig. H of model KM seems to increase with time 
indicating that this model is still evolving in time and 
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Fig. 6. The accretion rates of several quantities are plotted as a function of time for the highly supersonic (.Moo = 10) models JM 
and JS with a velocity gradient of 3%. The left panels contain the mass and angular momentum accretion rates, the right panels 
the specific angular momentum of the matter that is accreted. In the left panels, the straight horizontal lines show the analytical 
mass accretion rates: dotted is the Hoyle-Lyttleton rate (Eq. 1 in Ruffert 1994a), solid is the Bondi-Hoyle approximation formula 
(Eq. 3 in Ruffert 1994a) and half that value. The upper solid bold curve represents the numerically calculated mass accretion 
rate. The lower three curves of the left panels trace the x (dotted), y (thin solid) and z (bold solid) component of the angular 
momentum accretion rate. The same components apply to the right panels. The horizontal line in the right Panel of model JM 
shows the analytic specific angular momentum value as given by Eq. (pi). It is outside the range of the plot for model JS. 



Fig. 7. Contour plots showing snapshots of the density together with the flow pattern in a plane containing the center of the 
accretor for all models with a velocity gradient of 20%. The contour lines are spaced logarithmically in intervals of 0.1 dex. The 
bold contour levels are labeled with their respective values (0.01 or 1.0). The dashed contour delimits supersonic from subsonic 
regions. The time of the snapshot together with the velocity scale is given in the legend in the upper right hand corner of each 
panel. 
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Fig. 8. The accretion rates of several quantities are plotted as a function of time for the moderately supersonic (A^oo=3) 
models KM and KS with a velocity gradient of 20%. The left panels contain the mass and angular momentum accretion rates, 
the right panels the specific angular momentum of the matter that is accreted. In the left panels, the straight horizontal lines 
show the analytical mass accretion rates: dotted is the Hoyle-Lyttleton rate (Eq. 1 in Ruffert 1994a), solid is the Bondi-Hoyle 
approximation formula (Eq. 3 in Ruffert 1994a) and half that value. The upper solid bold curve represents the numerically 
calculated mass accretion rate. The lower three curves of the left panels trace the x (dotted), y (thin solid) and z (bold solid) 
component of the angular momentum accretion rate. The same components apply to the right panels. The value (-0.9) of the 
specific angular momentum as given by Eq. (pi) is outside the range of the plot for both models. 



has not yet reached a steady mean state. These pan- 
els should be compared to the analogous panels (left in 
Fig. U for the models IM and IS. One notices that both 
the mean and the amplitude of the fluctuations is larger 
in model KM than in model IM Although the mean mass 
accretion rate of model KS is lower than the mean in 
model IS this might be a transient effect, since toward 
the end of model IS the mass accretion rate is of the same 
magnitude as in model KS. Thus for the larger accretors 
(M-models) the vortex allows more mass to be accreted, 
while for the small accretors (S-models) the difference is 
negligible. Both models KM and KS have one difference in 



common compared to model IM and IS respectively: the 
unstable flow, which manifests itself e.g. via a fluctuat- 
ing mass accretion rate, begins much faster for the mod- 
els with a large gradient, models KM and KS. Also the 
large-scale fluctuations of the specific angular momentum 
appear at roughly t « 1 in models KM and KS, while it 
takes until is; 4 in models IM and IS. 

5.2. Highly supersonic accretors, Mach 10 

The accretion rates of several quantities for the highly su- 
personic models LM and LS are shown in Fig. ^. The flue- 
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Fig. 9. The accretion rates of several quantities are plotted as a function of time for the highly supersonic (.Moo = 10) models LM 
and LS with a velocity gradient of 20%. The left panels contain the mass and angular momentum accretion rates, the right panels 
the specific angular momentum of the matter that is accreted. In the left panels, the straight horizontal lines show the analytical 
mass accretion rates: dotted is the Hoyle-Lyttleton rate (Eq. 1 in Ruffert 1994a), solid is the Bondi-Hoyle approximation formula 
(Eq. 3 in Ruffert 1994a) and half that value. The upper solid bold curve represents the numerically calculated mass accretion 
rate. The lower three curves of the left panels trace the x (dotted), y (thin solid) and z (bold solid) component of the angular 
momentum accretion rate. The same components apply to the right panels. The value (-3.0) of the specific angular momentum 
as given by Eq. (^) is outside the range of the plot for both models. 



tuation amplitude of the mass accretion rate of model LM 
is larger than both the amplitudes of model JM (top left 
panel in Fig. ^J) and of model FM (top left panel of Fig. 9 
in Ruffert 1994), however the mean seems roughly equal. 
On the whole model LM looks more unstable and active 
than the models JM and FM with small or no gradients. 
The mass accretion rate of model LS does not seem to de- 
cline constantly during the first two time units as was the 
case for model JS. 

The trend that the model with the smaller accretor 
(in this case model LS compared to model LM) has the 
z-component shifted closer to the zero line is repeated 



also for the highly supersonic models with large gradients. 
The fluctuations around the mean of the z-componcnt has 
roughly the same amplitude and frequency as the fluctu- 
ations of the x and y-components around zero. 

6. Results of models with index 4/3 

As a check for the numerics two models were run with the 
same physical initial and boundary conditions, but differ- 
ing in resolution of the accretor: model SS used 9 grids, 
model ST used 10 grids (cf. Tablejllfor the other parame- 
ters) . The two left panels of Fig. [jlj show contour plots of 
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Fig. 10. Contour plots showing snapshots of the density to- 
gether with the flow pattern in a plane containing the center of 
the accretor for models SS, ST and RL. The contour lines are 
spaced logarithmically in intervals of 0.1 dex. The bold con- 
tour levels are labeled with their respective values. The dashed 
contour delimits supersonic from subsonic regions. The time 
of the snapshot together with the velocity scale is given in the 
legend in the upper right hand corner of each panel. 



the density distribution, while the acretion rates of mass 
and angular momentum can be found in Fig. [H]. Espe- 
cially the temporal evolution of the accretion rates con- 
firm that also for 7 = 4/3 using grids to 9 levels deep 
is sufficient: within the time in which the two models SS 
and ST overlap, the mass accretion rate and the angular 
momentum accretion agree to within a few percent. 

Models SS and ST (Fig. [ll]) are equivalent to mod- 
els IS and IT, (Fig. g) so a comparison should show what 
effects are due to the different adiabatic index 7. Addi- 
tionally, model SS can be compared to model CS shown 
in Fig. 6 of Ruffert (1995). Both the mean and the ampli- 
tude of the mass accretion rate of model SS are similar to 
the ones of model CS, and even the angular momentum 
accretion fluctuations do not show any striking difference. 
When comparing model SS to model IS one notices that 
the mass accretion rate proceeds in a much more unstable 
way in model SS. The rate of mass accretion as well as the 
specific angular momentum of the ^-component are larger 
in model SS. So much so, that in a short burst model ST 
actually reaches the analytically estimated value (middle 
right panel of Fig. [ll]) contrary to model IS or IT. Thus 
we conclude that the models with 7 = 4/3 accrete in a 
less stable way than their 7 = 5/3 equivalents, which is 
the opposite from what had been observed in the sim- 
ulations without gradients (Ruffert 1995 and references 
therein). More 7 = 4/3 models are necessary before any 
more systematic statements can be made. 

7. Results of model RL 

Model RL is extreme in two aspects: firstly the radius of 
the accretor is one accretion radius in size and thus very 
large, and secondly the flow speed upstream of the accretor 
is subsonic, Mach 0.6. Additionally the gradient was cho- 



"1 .2 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 



1.0 

CD 

-1— ' 

0.8 

o 

"-(— ' 

CD 
i__ 

o 
o 

CO 

I 0.4 

CO 



0.6 



0.2 
0.0 



Ma=3 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



-infty -2.0 -1.5 -1.0 -0.5 0.0 
lg gradient strength 



Fig. 12. Mass accretion rates (units: Mbh) are shown as a 
function of the strength of the velocity gradient: 20% and 3% 
are the results from this work, while the values for models with- 
out gradient (at the x-axis position "-infty") are taken from 
Ruffert (1994a) and Ruffert & Arnett (1994). Diamonds (o) 
denote models in which the accretor has a radius of 0.02 Ra, 
triangles (A) models with 0.1 Ra- The large bold symbols be- 
long to models with a speed of Moo = 10, while the smaller 
symbols belong to models with .Moc = 3. The accretor radius 
and Mach number are also written near each set of points. 
All models have 7 = 5/3, except for the "star" which denotes 
model SS. The error bars extending from the symbols indicate 
one standard deviation from the mean (S in Table [i]). Some 
points were slightly shifted horizontally to be able to discern 
the error bars. 
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Fig. 11. The accretion rates of several quantities are plotted as a function of time for models SS, ST and RL. The left panels 
contain the mass and angular momentum accretion rates, the right panels the specific angular momentum of the matter that is 
accreted. In the left panels, the straight horizontal lines show the analytical mass accretion rates: dotted is the Hoyle-Lyttleton 
rate (Eq. 1 in Ruffert 1994a), solid is the Bondi-Hoyle approximation formula (Eq. 3 in Ruffert 1994a) and half that value. The 
upper solid bold curve represents the numerically calculated mass accretion rate. The lower three curves of the left panels trace 
the x (dotted), y (thin solid) and z (bold solid) component of the angular momentum accretion rate. The same components 
apply to the right panels. The horizontal line in the right Panel of model ST shows the specific angular momentum value as 
given by Eq. (17]) . It is outside the range of the plot for models SS and RL. 
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Fig. 13. The relative mass fluctuations, i.e. the standard 
deviation S divided by the average mass accretion rate M 
(cf. Table |l|), is shown as a function of the strength of the ve- 
locity gradient: 20% and 3% are the results from this work, 
while the values for models without gradient (at the x-axis 
position "-infty") are taken from Ruffert (1994a) and Ruffert 

6 Arnett (1994). Diamonds (o) denote models in which the 
accretor has a radius of 0.02 Ra, triangles (A) models with 
0.1 Ra- The large bold symbols belong to models with a speed 
of Moo ~ 10, while the smaller symbols belong to models with 
Moo = 3. All models have 7 = 5/3, except for model SS (with 

7 = 4/3) denoted by a star (*). 



sen to be large (0.2). These conditions, although extreme, 
might, however, be present in stellar binary systems with 
a wind from a nova (K. Schenker, personal communica- 
tion), thus the model is not totally academic, but one has 
to take care in relating the physical size of the stars to the 
radius of the accretor used in the simulations. The con- 
tour plot showing the density distribution together with 
the instantaneous velocities is shown in Fig. [10)3, while 
the corresponding accretion rates are shown in the bot- 
tom panels of Fig. [ll]. The large scale of the contour plot 
shows the velocity gradient clearly. 

After t w 10 time units model RL has practically 
reached a stationary state in which mass is accreted at 
a constant rate (bottom left panel in Fig. [ll]. This be- 
haviour as well as the value of the mass accretion rate is 
very similar to model SL in Ruffert (1994b; top left panel 
of Fig. 2). The small fluctuations of the angular momen- 
tum visible in model SL are initial transients due to a 
small (3%) perturbation of the initial density distribution. 
I did not perturb the density in any model presented in 
this paper. This model RL is the only model presented in 
this work that reaches a constant state, which is probably 
due to the subsonic flow, sa has been reported elsewhere 
for models without gradients (e.g. Ruffert 1995). 

The main difference concerning accreted quantities be- 
tween model RL and all other models in this paper is 
the sign of the accreted angular momentum: it is positive 
(cf. Fig[ll]), while the z-component of the angular mo- 
mentum accreted of all other models is negative (see also 
Table |l|) . This is probably due to the fact that an accretor 
with such a large radius accretes mainly via its large geo- 
metric cross section, and not following the classic Bondi- 
Hoyle-Lyttleton scenario involving gravitational focussing, 
accretion mainly from the downstream side, etc. (the low 
Mach number might also play a role). Thus the matter 
is directly absorbed from the incoming stream without 
first being "processed" though a shock, taking with it the 
angular momentum it had upstream: this produces a anti- 
clockwise rotating vortex around the accretor (contrary to 
the clockwise vortex mentioned in Sect. ]|). 

8. Analysis of Results 

8.1. Mass accretion rate 



The mass accretion rates obtained in this work for all 
models except IT, RL, and ST are collected in Fig. [l2| to- 
gether with the amplitude of their fluctuations (one stan- 
dard deviation). This figure indicates that, to first order, 
the mass accretion rate is independent of e v , even for the 
large e v — 0.2. Although the mean mass accretion rate 
does vary slightly with Mach number and accretor radius, 
even in units of Mbh, the variation of the rates across 
different velocity gradients e v remains within the fluctu- 
ations of the unstable flow. This is only in partial agree- 
ment with Fig. |j Up to £ v = 0.1 no variation is expected 
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from Fig. [l| if all mass in the accretion cylinder is actu- 
ally accreted (which is an assumption that enters when 
deriving Eq. (||), etc.). However, since the mass accretion 
rate seems to remain unchanged even for the models with 
e v = 0.2 (contrary to Fig. Q), we conclude that for these 
large velocity gradients not all matter in the accretion 
cylinder is accreted any longer. 

The relative mass fluctuations, i.e. the standard devi- 
ation S divided by the average mass accretion rate M 
(cf. Table 0), have been collected in Fig. Although 
some models cluster around a relative fluctuation of 5%- 
10% while others are around 15%-25%, it is not clear 
which combination of parameters is responsable for this 
division. 

Two general statements can be made. When, starting 
from axisymmetric models, slightly increasing in the ve- 
locity gradient to a few percent, the relative accretion rate 
fluctuation either remains unchanged or decreases. In the 
axisymmetric case, any eddy will produce a fluctuation 
of the accretion rates. As long as the velocity gradient is 
small, the vortex generated by the the incoming angular 
momentum around the accretor is of the same strength 
as the eddies and it might be able to stabilize the flow 
around the accretor. When further increasing the velocity 
gradient the relative fluctuations increase strongly. So the 
stabilizing effect is lost indicating that the vortex itself 
contributes to the eddies and the fluctuations. 

8.2. Specific angular momentum 

Assuming a vortex flowing with Kepler velocity V just 
above the accretor's surface with radius R+, the specific 
angular momentum of such a vortex is 



10.00 



l s = R*V = ^R^/RaMocRac 



(10) 



The term y/2 has erroneously been omitted in my previ- 
ous works (e.g. Ruffert, 1996). Although for short periods 
of time the specific angular momentum can exceed l s , it 
is difficult to imagine how accreted matter can on average 
(temporal) exceed this value. This implies that smaller ob- 
jects (smaller R+) can accrete only smaller specific angular 
momenta, which goes to zero like \/R~l. 

In Fig. [14], I plot for several models the numerically ob- 
tained quantities l z along with the amplitude of the fluc- 
tuations (one standard deviation, <r z ), which can be found 
in Table 0. These are plotted in units of l s (Eq. (|l0|)). Ad- 
ditionally, above the diamonds denoting the above men- 
tioned ratio l z /l s , I plot, using plus-signs, the values that 
one expects for the analytically estimated quantity j z from 
Eq. ^ and denoting by squares, j z from the semi-numerical 
estimate, Eq. (||). 

Several trends can easily be noticed in Fig. E3. 
Model JM seems to be well below the general trend, in- 
dicating that the simulation was not evolved for long 
enough; I will not include this model in the following dis- 
cussion. The four "K" and "L" models form a fairly ho- 
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Fig. 14. The average specific angular momentum (units: l s , 
Kepler velocity vortex at surface of accretor, as given by 
Eq. (JlOT) ) is shown for most models by diamond symbols (l z in 
TableTO). The "error bars" extending from the symbols indi- 
cate one standard deviation from the mean (cr z in Table 0). The 
long error bars extending to the bottom axis are an indication 
that the fluctuations of the respective model are so large, that 
the specific angular momentum changes sign from time to time. 
The plus signs above the diamonds indicate the specific angu- 
lar momentum j z according to the Shapiro & Lightman (1976) 
prescription, Eq. (^), while the squares denote the values 
taken from the semi-numerical estimate, Eq. (^) and Fig. 
All models have 7 = 5/3, except for model SS (7 = 4/3). The 
star (*) at the position of I M i s the value taken from Ishii et 
al. (1993) (see text in Sect 



mogeneous group accreting roughly 0.3 of the Kepler spe- 
cific angular momentum. For the "I"- and "J" -models this 
fraction is roughly 0.1, thus confirming that for models 
with smaller gradients, the vortex around the accretor is 
less pronounced. These two groups vary less among them- 
selves than the variation one would expect if the analyti- 
cal estimates Eq. (0) (plus signs) or Eq. @ (squares) were 
valid. Thus, when estimating the specific angular momen- 
tum one should be guided by the Kepler- values Eq. (fio|). 
When applying Eq. @ one should bear in mind the al- 
lowable parameter range: small gradients, supersonic flow 
and small accretors. A good counter example is model RL, 
which exhibits a constant state and the sign of the accreted 
angular momentum is opposite to Eq. (^). 

The smaller lever arm acting in models with smaller 
accretors is included in Eq. jl0|). Still, the "S" -cases have a 
slightly smaller value of l z /l s compared to the "M"-cases. 
The reduction is, however, not uniform for all models; the 
models with large gradients show reductions of at most 
a factor of 2, while the other models have a factor of 3. 
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So small accretors impede high specific angular momen- 
tum accretion in some additional way than only via their 
smaller lever arm. At the distance of the surface (= ra- 
dius of the accretor), matter seems to move in eddies at 
some fraction of the Kepler-speeds. This fraction is depen- 
dent more on the velocity gradient e v than on the size of 
the accretor or the Mach number. Recall, that from Fig. [I] 
one would expect the models with large gradients ("K" 
and "L" ) to accrete specific angular momenta a factor of 
2 larger than the analogous models with small gradients 
("I" and "J"). This is in contradiction to what is shown in 
Fig. confirming again that the assumption of accreting 
everything from the accretion cylinder is not correct for 
large gradients. 

When the fluctuations of the specific angular momen- 
tum are larger than its mean, then the accreted angu- 
lar momentum can change sign, indicating a reversal of 
the rotation direction of a disk around the accretor. This 
will most easily be attained for models with a small gra- 
dient, since the fluctuations need not be large in these 
cases. These models have their "error bars" extending 
completely down to the x-axis in Fig. g For one of these 
models, model SS, Fig. [l5| shows the reversal of the disk 
surrounding the accretor. The "normal" rotation direc- 
tion is shown in Fig. |l5| c for model KS which has a large 
gradient e v = 0.2. From the bottom right panel in Fig. [s] 
one can see that the specific angular momentum never 
changes sign indicating a strong regular flow around the 
accretor. This is visible in Fig. [15|e: the disk (region of az- 
imuthal flow) extends to a distance of at least 0.3 Ra and 
is, however, slightly eccentric. In contrast, the disk around 
the models with small gradient, e.g. model SS, is smaller: 
roughly 0.1 Ra- Figs. |l5| a to d show this small disk al- 
ternating between clockwise and anti-clockwise rotation. 
Shocks appear when matter originating outside the disk 
is accreted in the opposite direction of the current disk 
rotation, e.g. shortly after the disk has been counterrotat- 
ing (in Fig. p"5|a), matter from downstream forces the disk 
to rotate in anti-clockwise direction (in Fig. |l5|b) . In this 
figure a shock is visible at (x w —0.05, y w —0.05): note 
the velocity discontinuity, the change from supersonic to 
subsonic (dashed line) and the increase in density (darker 
shades) when going with the flow in anti-clockwise direc- 
tion at (x « -0.05, y w -0.05). 

Shocks imply generation of entropy, thus if shocks ap- 
pear more often when the direction of rotation of the disk 
is reversed one would expect that matter with higher spe- 
cific entropy is accreted during phases when the specific 
angular momentum of the matter changes sign. To find 
such a possible correlation, I draw in Fig. [l6] a dot con- 
necting the two quantities for every second time step of 
the numerical simulation. This was done only for the four 
models JM, JS, IS, and SS which exhibited a change in 
sign of the specific angular momentum (see e.g. Fig [li]) . 
The two models JM and IS (Fig. [T^a and c) have an only 
small intrinsic scatter of the specific entropy accreted 



(roughly 10%) indicating that the fluctuations visible in 
Figs. [| and|6| do not generate equally fluctuating shock 
structures. For model JM the entropy does not at all 
seem to correlate with the angular momentum, while for 
model IS a marginal indication exists that the specific en- 
tropy is slightly higher for larger (=more positive) spe- 
cific angular momenta. On the other hand, the highest 
entropies appear at the most negative momenta. In con- 
trast, models JS and SS exhibit a relatively large scatter of 
the specific entropy (roughly 40% and 25%, respectively) 
and fairly clear correlations. In model JS (Fig. |l6|b) high 
entropy material is accreted preferentially when the an- 
gular momentum is small (around zero). This result con- 
firms what has been described at the beginning of this 
paragraph. The correlation in Model SS (Fig. [L6|d) is dif- 
ferent again: no material with low specific entropy (less 
than say 4.9) is accreted when the specific angular mo- 
mentum is positive. Thus the disk seems to be in con- 
stant turmoil (with many shocks) when the disk rotates 
in anti-clockwise direction (which is contrary to the "nor- 
mal" direction, cf. beginning of this Sect. 8.2). However, 
when the disk rotates in clockwise direction both high and 
low entropy material (more or less than 4.9) is present. 

In Fig. |l7] the rate at which mass is accreted is plotted 
versus the specific angular momentum. Analogously to the 
differing correlations of the entropy, the four models JM, 
JS, IS, and SS display different behaviours. The correla- 
tion in model SS is clearest: the mass accretion rate is 
highest when the specific angular momentum is around 
zero and the rate decreases for more positive and more 
negative values of the specific angular momentum. A sim- 
ilar, but less clear, trend can be discerned for model JM 
(Fig. [It]). Obviously, when the flow does not rotate around 
the accretor it falls down the potential to the surface of 
the accretor and can thus be absorbed more efficiently. 

The maximum rate at which mass is accreted also 
decreases with increasing magnitude of specific angular 
momentum (independent of sign) in models JS and IS 
(Figs. p~7| b and c). However, the smallest mass accretion 
rates are scattered fairly uniformly along all momenta. 

8.3. Comparison with previous works 

Two of the previously published fwo-dimensional simula- 
tions mentioned in the introductory chapter Sect. |l| in- 
vestigate a velocity gradient (as was done here); Taam & 
Fryxell (1989) and Anzer et al. (1987). The most impor- 
tant parameters of the first work are: the radius of the 
accretor is R= 0.037 Ra, the adiabatic index is 7= 4/3, 
the velocity gradient is e = 0.005 or e = 0.0625, with a 
Mach number of — A or Ai — 12. Thus a model most 
similar to these conditions is model SS (cf. Table 0). When 
comparing the top two panels of Fig. |ll| with the equiva- 
lent Figs. 10 and 11 in Taam & Fryxell (1989) one notices 
one main difference: model SS does not show the "flar- 
ing events" described by Taam & Fryxell (1989). These 



18 



M. Ruffert: Non-axisymmetric wind-accretion simulations 



Fig. 15. A closeup view of the matter and velocity distri- 
bution around the accretor for two models with small accretor 
(R = 0.02 -Ra). Arrows are overlayed to show the instantaneous 
velocity of matter. The contour lines are spaced logarithmically 
in intervals of 0.5 dex for model SS and 0.2 dex for model KS. 
The bold contour levels are sometimes labeled with their re- 
spective values (0.01 and 1.0). The dashed contour delimits 
supersonic from subsonic regions. The time of the snapshot 
together with the velocity scale is given in the legend in the 
upper right hand corner of each panel. 



flaring events are due to the collapse of otherwise fairly 
stable disks around the accretor. In quasi-regular inter- 
vals the disk changes its direction of rotation and dur- 
ing these inversions the mass of the disk is accreted. This 
yields short episodes of very high accretion rates, termed 
flaring events by Taam & Fryxell (1989). In model SS the 
disk is much less stable than in Sequence 2 of Taam & 
Fryxell (1989), consequently the buildup and collapse of 
the disk is much more eratic, so no flaring events of the 
same magnitude as in Taam & Fryxell (1989) is seen in 
model SS. That the disk is so stable in the simulation by 
Taam & Fryxell (1989) is due to the fact that their calcula- 
tion is two-dimensional, contrary to the three-dimensional 
models presented here. Once a disk is formed in a two- 
dimensional calculation hardly any matter can be accreted 
in radial direction. In three-dimensions, however, matter 
can still be accreted via the polar caps even if is disk is 
present in the equatorial plane. If this matter is focussed 
from above and below the disk, it can also act to disrupt 
the integrity of the disk, shortening the lifetime of the disk 
in three-dimensional calculations. 

The most important parameters of the Anzer et 
al. (1987) work are: the radius of the accretor is R — 
0.14 i?A, the adiabatic index is 7 = 1.5, the velocity gradi- 
ent is £ = 0.3 and the Mach number is M. = 3. Models IM 
or SS from the present work most closely resemble these 
parameters. In Sect. 4 of Anzer et al. (1987) they report 
finding a ratio of 0.22 between the numerically obtained 
(SPH simulation with 7500 particles) specific angular mo- 
mentum and the analytically estimated value. This value 
is within a factor of two of the results shown in Fig. [l4|. 

Of the five previously published f/iree-dimensional sim- 
ulations mentioned in the introduction, we will concen- 
trate our comparison to the results of Ishii et al. (1993). 
In the other cases the numerics is very coarse or question- 



able, e.g. few zones or particles, local time stepping not 
appropriate to the problem, etc. The parameters which 
were used in the "3D velocity inhomogeneous case" of Ishii 
et al. (1993) are the following. The radius of the accretor 
is R = 0.125 Ra, the adiabatic index is 7 = 1 (isother- 
mal), the velocity gradient is £ = 0.1 and their Mach num- 
ber is Al = 3. No model described in the present pa- 
per (cf . Table @) has such a low adiabatic index, but the 
other parameters are most closely covered by model IM 
or KM. The value that Ishii et al. (1993) obtain numeri- 
cally (24% is -0.11 devided by -0.45, these values are taken 
from Sect. 3.5 in Ishii et al. 1993) is represented by a star 
(*) in Fig. [lij. Their value (24%) must be used with cau- 
tion, since their model has an adiabatic index of 7 = 1 and 
their numerical resolution of the accretor is only two zones. 
Extrapolating from the models presented in this paper, I 
would expect a higher value, since the average specific an- 
gular momentum of model SS (smaller 7) is larger than of 
model IS (larger 7), and the models with larger accretors 
("M" -models) are larger still. 

9. Conclusions 

For the first time a comprehensive numerical three- 
dimensional study is presented of wind- accretion with a 
velocity gradient using a high resolution hydrodynamic 
code. I vary the following parameters: Mach number of 
the relative flow (Mach 3 and 10) , strength of the veloc- 
ity gradient perpendicular to this flow (3% and 20% over 
one accretion radius), radius of the accretor (0.02, 0.1 and 
1 accretion radius), and adiabatic index (5/3 and 4/3). 
The results are compared among the models with differ- 
ing parameters, to some previously published simulations, 
and also to the analytic estimates of the specific angular 
momentum of the matter that is accreted (Eq. (Q), which 
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Fig. 16. The specific entropy of accreted matter is plotted versus the specific angular momentum of this matter for models JM, 
JS, IS and SS. Each dot displays the two quantities at one moment in time. The times for which the dots are plotted are t>l, 
t\ J>2, ij;6, and i<;3, respectively. 



assumes that all angular momentum in the accretion cylin- 
der is actually accreted). 

1. All models with a small enough accrctor (with a size 
less or equal than 0.1 accretion radii) exhibit active 
unstable phases, very similar to the models without 
gradients. The accretion rates of mass, linear and an- 
gular momentum fluctuate with time, although not 
as strongly as published previously for 2D models 
(e.g. Fryxell & Taam 1988). Similarly to the 2D simu- 
lations, transient disks form around the accretor that 
alternate their direction of rotation with time. 

2. Depending on the model parameters, the average spe- 
cific angular momentum accreted is roughly between 
7% and 70% of the analytical estimate. For the models 
with small velocity gradients (3%) the accreted specific 
angular momentum is roughly a factor of 10 smaller 
than the value of a vortex with Kepler velocity around 



the surface of the accretor. This factor is roughly 3 for 
models with a large gradient of 20%. 

3. The mass accretion rates of all models with velocity 
gradients are equal, to within the fluctuation ampli- 
tudes, to the rates of the models without gradients 
(published previously) . 

4. The fluctuations of the mass accretion rate in the mod- 
els with small gradients (3%) are also similar to the val- 
ues of the models without gradients, while the models 
with large gradients (20%) exhibit larger fluctuations. 
So large gradients either amplify existing instability 
mechanisms or generate new ones. 

5. Marginal correlations are found, connecting the mass 
accretion rate, the specific angular momentum, and 
the specific entropy during the temporal evolution. 
The mass accretion rate is maximal when the spe- 
cific angular momentum is zero, while the specific en- 
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Fig. 17. The mass accretion rate is plotted versus the specific angular momentum of the accreted matter for models JM, JS, IS 
and SS. Each dot displays the two quantities at one moment in time. The times for which the dots are plotted are t > 1, t >2, 
tj>6, and £iJ3, respectively. 



tropy tends to be smaller when the disks are prograde 
(i.e. when the specific angular momentum is negative, 
in our units). 

Movies in mpeg format of the dynamical evolution 
of some models are available in the WWW at http: 
/ / www . mpa-gar ching . mpg . de/~mor/brila . html 
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